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Valley Polarization in Si (100) at Zero Magnetic Field 
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The valley splitting, which lifts the degeneracy of the lowest two valley states in a 
SiO2/Si(100)/SiO2 quantum well is examined through transport measurements. We demonstrate 
that the valley splitting can be observed directly as a step in the conductance defining a boundary 
between valley-unpolarized and polarized regions. This persists to well above liquid helium tem- 
perature and shows no dependence on magnetic field, indicating that single- particle valley splitting 
and valley-polarization exist in (100) silicon even at zero magnetic field. 

PACS numbers: 71.20.-b, 72.20.-1, 73,40.-c 
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Beyond its immense technological importance, sili- 
con continues to command considerable physical inter- 
est. This, on one hand, is driven by relentless technolog- 
ical advances producing ever- higher- quality and versatile 
structures such as Silicon on Insulator (SOI) allowing 
unexplored physical regimes to be examined, and com- 
plementary processing techniques facilitating new exper- 
imental schemes [ij. On the other hand, silicon possesses 
a number of unique properties, a very fundamental one 
of which is its electronic band structure. 

The conduction band dispersion relation of bulk sil- 
icon has 6 equivalent minima lying along three mutu- 
ally orthogonal axes of the Brilloin Zone. With two- 
dimensional confinement parallel to one of these axes, 
this six-fold degeneracy is lifted such that the two min- 
ima lying along the confinement axis become the lowest 
lying states due to anisotropy of effective mass 2j. Elec- 
trons in 2-dimensional structures such as Si(lOO) MOS- 
FETs (Metal Oxide Semiconductor Field Effect Tran- 
sistors) occupy these remaining two valleys such that 
each electron has both valley and spin degrees of free- 
dom. In recent years, this remaining valley degeneracy 
and valley splitting which lifts this degeneracy have ex- 
perienced a revival of interest, partly because valleys 
offer a relatively unexplored degree of freedom in rela- 
tion to metal-insulator transitions |3j and other many- 
body physi cs 14 I d I d Id, 1^ , but also due to new device 
possibilities|3,ll05lll| and as a source of decoherence in 
spin-based quantum computation |3, llO? II4 • 

Although effects indicating interactions between val- 
leys have been known for a long timejl^, they had been 
limited to structures involving vicinal planes which break 
the in-plane symmetry, or in (100) structures, special ex- 
perimental conditions of low temperature and magnetic 
field 0. Despite the long history, a consensus for the na- 
ture of this valley splitting has yet to be reached. In par- 
ticular, some studies I2I, Id, y^ll^ have pointed to a strong, 
almost linear, dependence on magnetic field, suggesting 
a many-body origin as opposed to a single-particle one, 
and its existence without magnetic field in Si(lOO) has 
not been experimentally confirmed. 



In this Letter, we show the first direct and vivid ev- 
idence for valley splitting of Si(lOO) in the absence of 
magnetic field. Recent high magnetic field experiments 
have found that the Si/Si02 interface in SIMOX (Sepa- 
ration by IMplantation of OXygen (a type of SOI)) struc- 
tures offer values of valley- splitting |l^ |l^ of at least a 
couple of meV, large compared with values found in con- 
ventional Si-MOS structures under equivalent conditions. 
Here, we show that the valley splitting in these structures 
can be increased to lO's of meV, and that this large valley 
splitting leads to a step feature in the conductance defin- 
ing a boundary between valley-unpolarized and valley- 
polarized regions. This allows us to verify the existence 
of valley-splitting in the absence of magnetic field. Fur- 
ther, we find its effects to be important well above liquid 
helium temperatures and that the valley splitting can be 
approximated by a linear dependence on electrical bias 
even at unprecedentedly large values found here. These 
findings strongly suggest its single-particle nature, and 
demonstrate its potential for extensive band structure 
control and for examining new physical regimes where 
valley splitting is comparable to, or greater than, other- 
wise dominant energy scales. 

The samples consist of SiO2/Si(100)/SiO2 quantum 
well transistors [Fig.l {inset)] fabricated on SIMOX sub- 
strates where the buried oxide is formed through ion- 
implantation followed by high temperature annealing 1^. 
Data presented here come from a nominally 8nm wide 
quantum well clad between thermal and buried oxide lay- 
ers 75nm and 380nm thick respectively, with large device 
dimensions of 30/im channel-width and 208/im between 
source and drain. The peak Hall mobility in these struc- 
tures is 0.8m^/Vs. Front and back-gate capacitances 
were found to be Cf = 485/iFm~^ and Cb = 92/iFm~^ 
respectively agreeing with sample dimensions. 

The source-drain current /sD was measured as a func- 
tion of front gate voltage Vyg at various values of back 
gate voltage I^g [Fig.l] (at a constant source volt- 
age (3mV) and constant temperature (4.2K), Gsd = 
/sD/3mV). With Vbg set at zero volts, Gsd shows be- 
havior characteristic of a standard silicon MOSFET at 



low temperature. At first, carriers induced by Vfq are 
localized and do not contribute to conduction. With the 
electron number density n increased, they begin to con- 
duct and GsD shows a sharp increase. When negative 
T^G is applied, the overall trace shifts to higher Ffg since 
the carrier concentration n also depends on Vbg accord- 
ing to n = Cf{Vfg - V^"") + Cb{Vbg - Vg^^) where V^'' 
and Vq^ are constant offsets. Also, since the electric 
field increases, a shallower slope in Gsd(^fg) results due 
to increased interface roughness scattering. Behavior at 
Vbg < is virtually identical to conventional MOSFETs 
with substrate bias 2]. 

With positive l^c, Gsd traces shift to lower T/pG and 
a minimum develops [A in Fig.l] due to the second sub- 
band, or alternatively worded, electron channels being 
present at both front and back interfaces [V^ . These two 
layers couple to form bonding and antibonding states, 
hereon referred to as "spatial subbands" , and the extent 
to which these should be described as two "electric" sub- 
bands or two spatially separated layers depends on the 
potential bias of the quantum well and self-consistent 
effects IQI |13- While interlayer coupling changes the 
wavefunction thereby affecting mobility jl^ |jj( , the ac- 
cessibility of the second subband allow further scattering 
processes 20, 21] which reduces the conductance. These 
inter-spatial-subband effects, which have already been 
examined in other systems [l^l, are expected to be in- 
volved in producing the minimum here. 

Positive l^G also shifts the wavefunction of the lower 
spatial subband closer to the buried-oxide interface with 
greater valley splitting 16]. With large positive values of 
yBG, the data reveal another feature [B in Fig.l]. We 
interpret this also to be caused by inter subband effects, 
but rather than spatial subbands, we argue that this fea- 
ture arises due to valley splitting leading to two valley 
subbands. At low Ffg, the conduction is solely due to 
the lower valley subband and the step feature represents 
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FIG. 1: (Color online) The source-drain conductance Gsd at 
4.2K as a function of Vfg at different values of Vbg- A and 
B mark features due to second spatial subband and upper 
valley subband occupation respectively. The inset shows the 
experimental setup with a TEM image of a SIMOX structure. 
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FIG. 2: Left: The conductance at Vbg = 60V, 4.2K. The 
lower two graphs show 1^* and 2^^ derivatives. The feature 
marked A is due to the second spatial subband. B marks 
the feature due to the upper valley while C marks the onset 
of conduction. Right: Gray-scale plot of the 2^^ derivative 
showing the evolution of A, B and C with (Vbg, Vfg)- White 
squares show calculated onsets of second spatial-subband oc- 
cupation [13. The abrupt horizontal line at Vfg = 2V is an 
experimental artifact. 



the onset of occupation of the upper valley [23. 

In order to examine the evolution of these features with 
Fbg, we doubly differentiate the GsdC^g) traces to ac- 
centuate the features [Fig. 2, left] and construct gray-scale 
plots as a function of (I^c'^g) [Fig. 2, right]. The ini- 
tial onset of conduction and the subsequent step both 
appear as peaks in the double-differential, and as white 
regions on the gray-scale plot. While peak C at low Vfq 
marks the onset of conduction, peak B defines the bound- 
ary between valley polarized and unpolarized regions. It 
is readily seen that with increased Fbg, these peaks sep- 
arate, and appear to continue separating in an approxi- 
mately linear manner. The valley-polarized region grows 
in range as V^g is increased indicating a linear increase in 
valley splitting A. This linear increase of A with V^g and 
hence electric field is qualitatively consistent with theo- 
retical predictions 2, 10] and previous magneto-transport 
estimates of the valley splitting 23]. 

To confirm that the structure arises due to valley- 
splitting, we performed measurements in magnetic field. 
Features similar to A,B and C at zero magnetic field 
can be seen at approximately the same (Fbg,^fg) at 
B = 5.5T (also marked A,B and C in Fig. 3(a)), de- 
marcating regions showing different patterns of Shub- 
nikov de Haas (SdH) oscillations. With positive V^g 
and small V^g (top-left of the graph) where the valley 
splitting A is known to be small, SdH oscillations can be 
seen with a period determined by the total filling factor 
^tot _ fiii/eB. Minima are seen whenever z/*°* is an in- 
teger multiple of 4, since 4 is the combined degeneracy 
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FIG. 3: (a) A^G^x^jAV^^ at 5.5T, 4.2K. (b) Same data as (a) 
but with different contrast. Solid squares joined by curves are 
estimated onsets of 2^^ spatial subband occupation |l7l \I\. 
Solid straight lines show calculated lines of even v^ and v~ 
while dotted lines represent zy*°* = M where i is an integer [23 • 
The region beneath z/~ = is completely valley-polarized. 
(c) Density of states at zero magnetic field with and without 
valley splitting A. (d) Resultant Landau levels. Thicker lines 
also correspond to valley subband edges. 



of valley and spin 2]. In contrast, in the region lying be- 
tween features B and C, the oscillations occur with half 
the periodicity of Az^*°* = 2, showing that the only re- 
maining degeneracy here is due to spin and the system is 
completely valley-polarized where the electrons only oc- 
cupy the lower valley subband. In the region above B, 
two sets of oscillations can be seen. One set describes 
lines rather parallel to B, verifying that B represents the 
onset of the upper valley subband [2J . 

To make estimates of A, we use a simple phenomeno- 
logical model. If each valley has identical and constant 
in-plane effective mass tti*, the density of states for each 
subband is I) = vri' vn^) j i\T? . When both valleys have 
non-zero occupation, the number of electrons in each 
valley n^ is then given by n^ = ^(^ ± DlS) in which 
superscripts + and - indicate lower and upper valley 
states respectively [Fig. 3(c)]. Since the directions along 
which the oscillations are occurring appear rather linear 
in (^FG, ^bg), we suppose a simple empirical approxima- 
tion: A = OL^n when ^n > where (5n = ne — np, n^ 
and np being electron densities contributed by the back 
and front gates, respectively. (We set A = for (5n < 
since it is negligibly small.) It then follows that 
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lower and upper valleys v^ and v~ respectively under 
magnetic field [Fig. 3(d)]. Solid lines in Fig. 3 (b) show 
regions where v^ and v~ are even, where the fitting pa- 
rameter has been set to a = 0.46meV/10^^m~^ chosen 
to fit the minimum along v~ = 2. Although this value 
of OL works reasonably well between i/~ = and v~ = 4, 
smaller values of around a = 0.40meV/10^^m~^ work 
better for higher filling factors between u~ = 8 and 10. 
This is likely to be due to the larger concentration lead- 
ing to more screening, reducing the effect of the external 
bias. In turn, it also suggests that greater values are 
present at lower concentrations. 

We have also performed measurements where the mag- 
netic field was swept keeping I^g constant at GOV which 
have shown that the SdH oscillations can be fitted to 
similar accuracy as the constant field data [Fig. 3] with- 
out any magnetic field dependence of A. This confirms 
that the valley splitting here is independent of magnetic 
field, corroborating the correspondence between the data 
at zero field [Fig. 2] and higher field [Fig. 3]. 

Fig. 4 shows GsD and its second derivative taken at 
zero magnetic field and I^g = 60V where we estimate 
A = 23meV (A/Zc^ = 270K) at the onset of upper val- 
ley occupation. The valley subband edges expected from 
the fit are also marked, suggesting that the actual bottom 
of the upper valley subband (n~ = 0) occurs at a lower 
value of VpG than the center of the peak in d^Gso/dFpQ. 
This is consistent with the interpretation that the step 
feature in Gsd arises as a consequence of suppression in 
mobility when the upper valley becomes energetically ac- 
cessible. The step-like increase occurs when this suppres- 
sion becomes weaker, after the onset of occupation, as 
the Fermi energy moves further above the subband edge. 
However, further work to establish the roles played by 
inter-valley interactions and the differences in mobility 
of the two individual valley subbands are required before 
a quantitative understanding of Gsd can be reached. 

Fig. 4 also shows that the conduction band edge ex- 
pected from the fit lies at lower I^g than the onset of 
conduction. In Fig. 3, it is also seen, that although 
lines of constant u~^ describe straight diagonal lines, the 
peak associated with the onset of conduction (C) curves 
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FIG. 4: {Color online) Temperature dependence at B = OT, 
Vbg = 60V. Left: Gsd- Right: d^GsD/d^FG- Arrows mark 
valley subband edges extracted using a — 0.46meV/10^^m~^. 



toward higher Vpc with increased I^bg- This can be at- 
tributed to increase in the number of localized states with 
Vbg a.s the disorder is increased, due in turn to the elec- 
trons being pressed against the Si-Si02 interface. With 
increased temperature, the peak in the second differential 
shifts toward the predicted conduction band edge. 

At higher temperatures, the feature associated with 
the upper valley can still be discerned directly in Gsd at 
33K, but it is not so clear at 77K. However, the double 
differential does show the second peak even at 77K. At 
such a temperature and above, the conduction globally 
becomes strongly influenced by phonon scattering and is 
no longer a good probe for detecting subband edges. 

We now discuss our findings in relation to recent ex- 
periments in other systems. Studies on conventional Si 
M0S|3 and Si/SiGe 8] have both shown very strong de- 
pendence of A on magnetic field. In Si/SiGe, the absolute 
values of A are extremely small (lO's of /ieV) suggesting 
that their phenomena belong to a completely different 
regime at very much smaller energy scales. In conven- 
tional Si MOS structures, A has relatively larger val- 
ues between 0.2 and 2meV, but these are still smaller 
than values addressed in this work. Furthermore, the 
reported field dependence comes from measurements of 
energy gaps in the quantized Hall regime 6, 14] and their 
results are likely to be heavily influenced by many-body 
effects dominant under their conditions. 

The results here suggest, therefore, that there exists a 
bare valley splitting arising from single particle interac- 
tions [lO|, independent of magnetic field. When it is of 
the order of a few meV and greater, it dwarfs many body 
contributions and there is very little magnetic field de- 
pendence. When the bare splitting is small, on the other 
hand, interaction effects become dominant and heavily 
influences the energy-gaps and its dependence on mag- 
netic field. 

There is also another important conclusion we can 
draw from the present results. In our previous experi- 
ments, we could only conclude that the valley splitting 
at the Si-buried oxide interface was much larger com- 
pared with the standard thermal-oxide interface at the 
front showing usual small values [l^. The linear depen- 
dence on electrical bias we find here, demonstrates that 
this rule is applicable even when the absolute value of 
A becomes orders of magnitude greater than those nor- 
mally observed. What changes from one type of interface 
to another is the prefactor a. 
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